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Abstract — We consider cooperative spectrum sensing for cog- 
nitive radios. We develop an energy efficient detector with low 
detection delay using sequential hypothesis testing. Sequential 
Probability Ratio Test (SPRT) is used at both the local nodes 
and the fusion center. We also analyse the performance of 
this algorithm and compare with the simulations. Modelling 
uncertainties in the distribution parameters are considered. Slow 
fading with and without perfect channel state information at the 
cognitive radios is taken into account. 

Keywords- Cognitive Radio, Spectrum Sensing, Cooperative 
Distributed Algorithm, SPRT. 

I. Introduction 

Cognitive Radio has evolved as a working solution for 
the scarcity of spectrum due to the proliferation of wire- 
less services. Cognitive Radios (CRs) access the spectrum 
licensed to other service providers opportunistically without 
interference to the existing communication services. For this 
the Cognitive users sense the spectrum to detect the usage 
of the channel by the primary (licensed) users. However due 
to the inherent transmission impairments of wireless channels 
and strict spectrum sensing requirements for Cognitive Radios 
ifTTl spectrum sensing has become one of the main challenges 
faced by them. 

Cooperative spectrum sensing ( 1201 . 11231 ) in which different 
cognitive radios interact with each other, is proposed as an 
answer to the problems caused by multipath fading, shadowing 
and hidden node problem in single node spectrum sensing 
methods. Also it improves the probability of false alarm and 
the probability of miss-detection. These are achieved via the 
exploitation of spatial diversity among the Cognitive users. 

Cooperative spectrum sensing can be either centralized or 
distributed ll23l . In the centralized algorithm a central unit 
gathers sensing data from the Cognitive Radios and identifies 
the spectrum usage ( l23l . lfT31 ). On the other hand, in the 
distributed case each secondary user collects observations, 
makes a local decision and sends to a fusion node to make the 
final decision. The information that is exchanged between the 
secondary users and the fusion node can be a soft decision 
(summary statistic) or a hard decision fl31 . Soft decisions 
can give better gains at the fusion center but also consume 
higher bandwidth at the control channels (used for sharing 
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information among secondary users). However hard decisions 
provide as good a performance as soft decisions when the 
number of cooperative users increases I0 . 

Spectrum sensing algorithms used at a node can use a 
fixed sample size (one shot) or sequential detection (Q, iflOll . 
|16|, l23l ). In case of fixed sample size detectors with the 
complete knowledge of primary signal, matched filter is the 
optimal detector 1 8 1 that maximises the SNR. When the only 
known apriori information is the noise power, the energy 
detector is optimal [8|. Sequential detection can provide better 
performance lfl2l . In the sequential approach one can consider 
detecting when a primary turns ON (or OFF) (change detec- 
tion) or just the hypothesis testing whether the primary is ON 
or OFF. Sequential change detection is well studied in (|2), 
iflOl . Ifl2l ). In sequential hypothesis testing (@, 0, [16)) one 
considers the case where the status of the primary channel is 
known to change very slowly, e.g., detecting occupancy of a 
TV transmission. Usage of idle TV bands by the Cognitive 
network is being targeted as the first application for cognitive 
radio. In this setup Walds' Sequential Probability Ratio Test 
(SPRT) provides the optimal performance for a single node 
( fl4l . Il22ll ). Il23l has an extensive survey of spectrum sensing 
methods. Other spectrum sensing schemes include methods 
based on higher order statistics fl3l . wavelet transforms lfl8ll 
and compressed sensing |fl9l . 

We use the sequential hypothesis testing framework in the 
cooperative setup. We use SPRT at each local node and again 
at the fusion center. This has been motivated by our previous 
algorithm, DualCUSUM used for distributed change detection. 
Thus we will call this algorithm DualSPRT. However this 
has been studied in ([9| and ifTBI ) as well. But unlike f|9l, 
1 16 1) we also provide theoretical analysis of this algorithm 
and consider the effect of fading in the channel between the 
primary and secondary nodes. We also model the receiver 
noise at the fusion node and use physical layer fusion to reduce 
the transmission time of the decisions by the local nodes to 
the fusion node. 

This paper is organised as follows. Section II describes 
the model. Section III starts with the DualSPRT algorithm. 
Simulation results and analysis are also provided in Section 
III. Then we consider the case where the SNRs are different 
at different Cognitive Radios. The received SNR may or may 



not be known to the CR nodes. In Section IV we introduce 
fading at the channel between the primary transmitter and the 
Cognitive Radios. The channel gains may not be available to 
the local secondary nodes. Section V concludes the paper. 

II. System Model 

We consider a Cognitive Radio system with one primary 
transmitter and L secondary users. The L nodes sense the 
channel to detect the spectral holes. The decisions made by the 
secondary users are transmitted to a fusion node via a Multiple 
Access Channel (MAC) for it to make a final decision. 

Let Xk,i be the observation made at secondary user I at 
time k. The {Xk,i, k > 1} are independent and identically 
distributed (iid). It is assumed that the observations are inde- 
pendent across Cognitive Radios. Based on {X n j, n < k} 
the secondary user I transmits Y^j to the fusion node. It is 
assumed that the secondary nodes are synchronised so that 
the fusion node receives Yj. = YliLi ^fc,i + ^fc> wnere {Zk} 
is iid receiver noise. The fusion center uses {Yk} and makes 
a decision. The observations {Xkj} depend on whether the 
primary is transmitting (Hypothesis Hi) or not (Hypothesis 
H ) as 

x _ ( Z k ,i, k = 1,2,..., under H Q 

k ' 1 \ hiSk + Zk,i, k = l,2, under Hi 

where hi is the channel gain of the I th user, Sk is the primary 
signal and Zk.i is the observation noise at the I th user at time 
k. We assume {Zk,i,k > 1} are iid. Let N be the time to 
decide on the hypothesis by the fusion node. We assume that 
N is much less than the coherence time of the channel so 
that the slow fading assumption is valid. This means that hi 
is random but remains constant during the spectrum sensing 
duration. 

The general problem is to develop a distributed algorithm 
in the above setup which solves the problem: 

xmnE DD =E[N\Hi] , (2) 

subject to Pfa < ot 

where Hi is the true hypothesis, i — {0,1} and Pfa is the 
probability of false alarm, i.e., probability of making a wrong 
decision. We will separately consider E[N\H{\ and E[N\Ho\- 
It is well known that for a single node case (L = 1) Wald's 
SPRT performs optimally in terms of reducing E[N\Hi] and 
E[N\Hq] for a given Pfa- Motivated by the good performance 
of DualCUSUM in (QJ, Q) and the optimality of SPRT for a 
single node, we propose using DualSPRT in the next section 
and study its performance. 

III. DualSPRT algorithm 

To explain the setup and analysis we start with the simple 
case, where the channel gains, hi=\ for all Vs. We will 
consider fading in the next section. DualSPRT is as follows: 

1) Secondary node, I, runs SPRT algorithm, 

W ,i = 

W k ,i = Wk-i,i + log [f h i (X k ,i) /h,i (X k ,i )] , k > 1 (3) 



where fu is the density of Xk,i under Hi and/o.; is the 
density of Xk,i under Ho. 

2) Secondary node I transmits a constant bi at time k if 
W k ,i > 71 or transmits b when W k ,i < Jo, i.e., Y k .i = 

&ll{Wk,i>7i} + & ol{W fc , i < 7o } 

where 70 < < 71 and 1a denotes the indicator 
function of set A. Parameters 61,60,71,70 are chosen 
appropriately. 

3) Physical layer fusion is used at the Fusion Centre, i.e. , 
Yk = Ylt=i Y~k,i + Zk, where Zk is the iid noise at the 
fusion node. 

4) Finally, Fusion center runs SPRT: 

F k = F k -i+ log [gi(Y k )/ g (Y k )}, F =0, (4) 

where g is the density of Zk + /io, the MAC noise at 
the fusion node, and gi is the density of Zk + H\, /io 
and (tti being design parameters. 

5) The fusion center decides about the hypothesis at time 
N where 

N = inf{fc : F k > (ii or F k < /?o} 

and (3q < < /3i. The decision at time N is Hi if 
Fn > ft.) otherwise Ho- 
In order to have equal Pfa under both hypothesis, we choose 

7i = -7o = 7 and ft = -ft = ft 

Of course Pfa can be taken different under Ho or Hi 
by appropriately choosing 71, 70, ft, ft. Any prior in- 
formation available about Ho or Hi can be used to de- 
cide constants. Performance of this algorithm depends on 
(7i)70)ft)ft)^i)&0!A*i)Mo)- Also we choose these parame- 
ters such that the probability of false alarm, Pf a at local nodes 
is much lower than Pfa- A good set of parameters for given 
SNR values can be obtained from known results of SPRT. 

Deciding at local nodes and transmitting them to the fusion 
node reduces the transmission rate and transmit energy used 
by the local nodes in communication with the fusion node. 
Also, physical layer fusion in Step 3 reduces transmission 
time, but requires synchronisation of different local nodes. If 
synchronisation is not possible, then some other algorithm, 
e.g., TDMA can be used. 

DualSPRT (without physical layer synchronization and fu- 
sion receiver noise) has been shown to perform well in (|9;j, 
[16 1). In the rest of the following we analyse the performance 
under our setup. 

A. Performance Analysis 

We first provide the analysis for Ebd and then for Pp^.The 
analysis for Edd is similar to that of DualCUSUM in Q. 
For simplicity, in the following we will take 71 = —70 = 7, 
ft = — ft = ft Mi = ^ Mo = M an d ft = —bo = 1. Then 
Pfa under the two hypothesis is same. 



Eod Analysis: At the fusion node F k crosses /3 under 
Hi when a sufficient number of local nodes transmit b\. 
The dominant event occurs when the number of local nodes 
transmitting are such that the mean drift of the random walk 
Fk will just have turned positive. In the following we find the 
mean time to this event and then the time to cross /3 after 
this. The Eod analysis is same under hypothesis Ho and H\. 
Hence we provide the analysis for H\. 

At secondary node I SPRT {Wk,i,k > 0} is a ran- 
dom walk. Let d = E Hl [log (A //„ (X M ))] , <J 2 = 
Var[\og [/i {X k>l ) //„ (X kjl )]]. We know S > 0. The time r 7 
for W k at each local node to cross the threshold 7 satisfies 
E[t^\ ~ 7/(5 for large values of 7 (needed for small Pfa- 
Then by central limit theorem we can show that at each node 



P Hl (FA before h) 



Now, as in J7J, we can show that, 



Eod ~ E[tj 



(3-F, 



(5) 



(6) 



where Sj is the drift of the fusion center SPRT, Fk when j 
local nodes are transmitting, tj is the point at which the drift 
of Fk changes from Sj-i to Sj, Fj — E[F t -i], the mean 
value of Fk just before transition epoch tj, and 



j = min{i : Si > and 



(3-F 



< E[t i+1 ] - E[U}}. 



An iterative method is proposed [2| to calculate E[ti] and Fj 
in an efficient manner. 

For the above analysis for Edd we followed the analysis 
of DualCUSUM in [7 |. However there are some difference in 
the SPRT at the fusion center here from the DualCUSUM in 
0. But comparison with simulations show that we will get 
an acceptable approximations. 

Pfa Analysis: It can be easily verified that tk, defined 
earlier is the k th order statistics of L iid random variables, 
77./ (first passage time to threshold 7 by the I th node, whose 
probability density function is given in (0) . Then Pfa when 
Hi is the true hypothesis is given by, 

PH t {False alarm) = PH t {False alarm before ti) (7) 

+PH 1 {False alarm between tiand £2) 

+ Pn 1 (False alarm between ti and £3) + .... 

One expects that the first term in (O should be the 
dominant term. This is because Pj a is much smaller than 
Pfa and hence after t\, the drift of Fk will be more positive. 
Therefore the probability of false alarm goes down. We have 
verified this from simulations also. Hence we focus on the 
first term. 

Let S k = log [ 9x (Y k ) /g (Y k )] and 9 = p/2u. 

Therefore F k = Si + S2 + ■ ■ ■ + Sk ■ Every Si , 1 < i < k has 
a common term 2u (in case of Gaussian gi and go) , thus 
changing the threshold to 9 = (3/2u. Then 
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p[s k < -cifp^.i-e + c}dc 



(l-2P[F k -i < -6>]Vl-$ tl (fe) 



where $ fl is the Cumulative Distribution Function of t\. 
As we are considering only {Fk,k < ti}, we remove the 
dependencies on t\. (A) is because of the Markov property of 
the random walk. (B) is due to the inequality, 

P[supF fc > 9] < 2P[F n > 9} 

for the Gaussian random walk F k [4|. Similary we can write 
an upper bound by replacing P[n^Z^{F„ > —9}] with 
P[Fk-i > —9]. In Table U we compare the lower bound 
on Pfa with the simulation results. We can make this lower 
bound tighter if we do the same set of analysis for the Gaussian 
random walk between ti and ti with appropriate changes and 
add to the results we already obtained. 

B. Example 

We apply the DualSPRT on the following example and 
compare the Eoo and Pfa via analysis provided above with 
the simulation results. We assume that the pre-change distribu- 
tion fo and the post-change distribution fi are Gaussian with 
different means. This model is relevant when the noise and 
interference are log-normally distributed l20l . This is a useful 
model when Xk.i is the sum of energy of a large number of 
observations at the secondary node at low SNR. 

Parameters used for simulation are as follows: There are 5 
secondary nodes, (L = 5), f ~ 7V(0, 1) and fi - J\f(l, 1), 
where Af(a, b) denote Gaussian distribution with mean a and 
variance b. Also fo = foj and fx = fu for 1 < I < L, 71 = 

-7o = 7. Pi = ~/3o = P, Ml = -Mo = A 1 and h = -b = 1. 
The Pfa and the corresponding Eod are provided in Table 
U The parameters are chosen to provide good performance for 
the given Pfa- The table also provides the results obtained 
via analysis. 

C. Analysis for different SNRs 

The above analysis is for the case when Xk,i have the 
same distribution for different / under the hypothesis Ho and 



hyp 


PpASim. 


Pp^Anal. 


E DD Sim. 


E u £, Anal . 


HI 


0.00125 


0.0012 


15.6716 


16.4216 


HI 


0.01610 


0.0129 


13.928 


12.6913 


HO 


0.0613 


0.0497 


11.803 


10.583 


HO 


0.0031 


0.0027 


15.1766 


14.830 



TABLE I 

DualSPRT: Comparison of E dd and P fa obtained via analysis 
(lower bound on the dominating term) and simulation 

H\. However in practice the Xk,i for different local nodes I 
will often be different because their receiver noise can have 
different variances and / or the path losses from the primary 
transmitter to the secondary nodes can be different. The above 
analysis for this case needs slight changes for Edd as well 
as P FA . 

For the analysis of Edd one difference is that r 7 ./, I = 
1, . . . , L are no longer iid. Now the iterative scheme used in 
Section III A to calculate E t - and Fj does not work. Thus, 
knowing the minimum number of local nodes needed to make 
the mean drift of Fk positive (say it is i*), we compute the 
mean of the i* order statistics of the independent random 
variable r 7i ;, I = 1, . . . , L via J3]. Then we approximate the 
Edd by 

p_fE[t i . ] -Eit i ._ l] \ 
E[UA + " 7~~~~ ■ ( 8 > 



For Pfa analysis we need the distribution of the first order 
statistics t\ for T lt i,l — l...,L and then use the method 
proposed in Section III A. 

We provide an example to verify the accuracy of the 
performance analysis provided above. 

D. Example 

There are five secondary nodes with primary to secondary 
channel gain being 0, -1.5, -2.5, -4 and -6 dB respectively 
(corresponding post change means are 1, 0.84, 0.75, 0.63, 0.5). 
f ~ 7V(0, l),/ = fo.i for 1 < I < L. Table M provides the 
Edd and Pfa via analysis and simulations. We see a good 
match. 



PpASim. 


PpAAnal. 


EooSiva. 


Edd Anal. 


26.68e - 4 


27.51e-4 


36.028 


34.634 


18.78e - 4 


19.85e - 4 


44.319 


43.290 


36.30e - 4 


35.16e - 4 


27.770 


25.977 



TABLE II 

DualSPRT for different SNR' s between the primary and the 

SECONDARY USERS: COMPARISON OF Ertrt AND P F A OBTAINED VIA 
ANALYSIS AND SIMULATION. 

E. Different and unknown SNRs 

Next we consider the case where the received signal power 
is fixed but not known to the local Cognitive Radio nodes. 
This can happen if the transmit power of the primary is not 
known and / or there is unknown shadowing. Now we limit 
ourselves to the energy detector where the observations X^ : i 
are a summation of energy of N samples received by the I th 
Cognitive Radio node. Then for somewhat large N, the pre 
and post change distributions of Xk.i can be approximated 
by Gaussian distributions: fo.i ~ Af(o~ 2 , 2a A /N) and fu ~ 
M(Pi + a 2 ,2(Pi + a 2 ) 2 /N), where Pi is the received power 



at the I th CR node and noise Z k ,i ~ Af(0,a 2 ). Under low 
SNR conditions (Pi + a 2 ) « <r 4 and hence X^ i are Gaussian 
distributed with mean change under Ho and H\. Now taking 
Xk,i — v 2 as the data for the detection algorithm at the I th 
node, since Pi is unknown we can formulate this problem as 
a sequential hypothesis testing problem with 

H Q :9 = 0;H 1 :9>9i. (9) 

where 9 is Pi and 9\ is appropriately chosen. 
The problem 

H : 9 < 9 ; H x : 9 > 6 1 , (10) 
subject to the error constraints 

Pe{rejectH a } < a for 9 < 9 (11) 

PelrejectH^ < j3 for 9 > 9 X 

for exponential family of distributions is well studied in ( ifTTIl . 
lfT2ll ). The following algorithm of Lai is asymptotically Bayes 
optimal ifTTI and hence we use it at the local nodes instead of 
SPRT. Let 9 6 A = [ai,o 2 ]. Define 



max 



log 



feAXk) 

N(g,c) = inf {n : W n .i > g(nc)} , 



(12) 



(13) 



where g() is a time varying threshold. Its approximate expres- 
sion is given in IfTTIl . At time N(g,c) decide upon Ho or Hi 
according as 



a N{g,c) < V Or °N{g.c) 

where 9* is obtained by solving 1(9* ,9 ) — 1(9* ,6{), and 
1(9, A) is the Kullback-Leibler information number. Also for 
Gaussian f and fi, 9 n = mos{ai, ram[S n /n, 02]}. 

The choice of 9\ in (0 affects the performance of E[N\Ho] 
and E[N\Hi] for the algorithm (|T2)-(fT3), where N = N(g, c). 
For our case where Ho : 9 = 0, unlike in (TTOb where Hq : 
9 < 0, E[N\Hq] largely depends upon the value 9%. As 9\ 
increases, E[N\Ho] decreases and E[N\Hi] increases. If Pi 6 
[P,P] for all I then a good choice of 9 U is (P - P)/2. 

In the distributed setup with received power at the local 
nodes unknown, the local nodes will use the Lai's algorithm 
mentioned above while the fusion node runs the SPRT. All 
other details remain same. We call this algorithm GLR-SPRT. 

The performance of GLR-SPRT is compared with Dual- 
SPRT (where the received powers are assumed known at the 
local nodes) for Example III D in Table |IjTJ Interestingly 
E[N\Hi] for GLR-SPRT is actually lower than for DualSPRT 
, but E[N\H ] is higher. 



> 



hyp 


Edd 


Pfa = 0.1 


Pfa = 0.05 


Pfa = 0.01 


HI 


DualSPRT 


2.06 


3.177 


5.264 


HI 


GLRSPRT 


1.425 


2.522 


4.857 


HO 


DualSPRT 


1.921 


3.074 


5.184 


HO 


GLRSPRT 


2.745 


3.852 


6.115 



TABLE III 

Comparison between GLRSPRT and DualSPRT for different 

SNR'S BETWEEN THE PRIMARY AND THE SECONDARY USERS. 



IV. Channel with Fading 

In this section we consider the system where the channels 
from the primary transmitter to the secondary nodes have 
fading (hi ^ 1). We assume slow fading, i.e., the channel 
coherence time is longer than the hypothesis testing time. We 
consider two cases, Case 1: the fading gain is known to the 
CR nodes. Case 2: the fading gain is not known to the CR 
nodes. 

When the fading gain hi is known to the I th secondary 
node then this case can be considered as the different SNR 
case studied in Section III C. Thus we only consider Case 2 
where the channel gain hi is not known to the I th node. 

We consider the energy detector setup of Section III E. 
However, Pi, the received signal power at the local node 
I is random. If the fading is Rayleigh distributed then Pi 
has exponential distribution. The hypothesis testing problem 
becomes 

H : /ci-AT^a 2 );^ : f u ~ Af(0, a 1 ) (14) 

where 9 is random with exponential distribution and a 2 is the 
variance of noise. We are not aware of this problem being 
handled via sequential hypothesis testing. However we use 
Lai's algorithm in Section III E where we take 9\ to be the 
median of the distribution of 6, such that P(9 > 9\) = 1/2. 
This seems a good choice for 9\ to compromise between 
E[N\H ] and E[N\Hi]. 

We use this algorithm on an example where a 2 = 1,9 = 
exp(l), Var(Zfe) = 1, and L = 5. The performance of this 
algorithm is compared with that of DualSPRT (with perfect 
channel state information) in Table [IV] (under Ho) and Table fVl 
(under Hi). The Edd and Pfa were computed by simulations 
each case by 100000 times and taking the average. We observe 
that under H\, for high Pfa this algorithm works better than 
DualSPRT with channel state information, but as Pfa de- 
creases DualSPRT becomes better and the difference increases. 
For Ho, GLRSPRT is always worse and the difference is 
almost constant. 



Edd 


Pfa = 0.1 


Pfa = 0.05 


P FA = 0.01 


DualSPRT 


1.669 


2.497 


4.753 


GLRSPRT 


3.191 


4.418 


7.294 



TABLE IV 

Comparison between GLRSPRT and DualSPRT with slow-fading 

BETWEEN PRIMARY AND SECONDARY USER UNDER HO. ENERGY 
DETECTION STATISTIC IS USED AT THE SECONDARY NODES 



Edd 


Pfa = 0.1 


Pfa = 0.08 


Pfa = 0.06 


DualSPRT 


1.74 


1.854 


2.417 


GLRSPRT 


1.62 


3.065 


5.42 



TABLE V 

Comparison between GLRSPRT and DualSPRT with slow-fading 

BETWEEN PRIMARY AND SECONDARY USER UNDER H 1 .ENERGY 

DETECTION statistic is used at the secondary nodes 
V. CONCLUSIONS AND FUTURE WORK 

We have proposed an energy efficient, distributed coopera- 
tive spectrum sensing technique, DualSPRT which uses SPRT 
at the cognitive radios as well as at the fusion center. We also 
provide analysis of DualSPRT. Next we modify the algorithm 



so as to be able to detect when the received SNR is not known 
and when there is slow fading channels between the primary 
and the secondary nodes. Future work should consider analysis 
of the GLR algorithms and optimising over the current setup. 
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